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ABSTRACT 


Akima  cubic  spline  and  spatial  Fourier  series  (SFS)  techniques  for  modeling  ocean 
bottom  contours  from  bathymetric  data  were  comparatively  analyzed.  SFS  methods  en¬ 
countered  difficulty  in  ocean  bottom  reconstruction  despite  several  enhancements  to  the 
fundamental  technique.  Both  Akima  cubic  spline  and  SFS  approaches  proved  unsatis¬ 
factory'  in  reproducing  first  and  second-order  derivatives  for  several  arbitrarily  shaped 
ocean  bottom  contours. 

A  simple  reflection  angle  algorithm  for  arbitrary  one-dimensional  bottom  models 
was  examined  and  found  to  be  accurate.  A  graphical  demonstration  of  acoustic  ray 
interaction  with  the  ocean  bottom  using  a  variety  of  mathematical  ocean  bottom  models 
and  the  reflection  angle  algorithm  was  conducted. 
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I.  INTRODUCTION 


The  primary  endeavour  of  this  thesis  research  was  the  comparative  evaluation  of 
existing  numerical  cubic  spline  techniques  with  spatial  Fourier  series  (SFS)  methods  for 
the  interpolative  modeling  of  real,  discrete,  ocean  bathymetric  data.  The  Akima  cubic 
spline  was  chosen  for  its  reported  ability  to  produce  ocean  bottom  contour  curves  which 
pass  through  the  original  data  while  minimizing  artificial  oscillations  in  the  interpolation. 
The  Akima  cubic  spline  logic  was  implemented  via  the  IMSL  Version  10  double  preci¬ 
sion  computer  program  DCSAKM.  First  and  second-order  derivatives  of  the  resulting 
spline  were  calculated  from  a  second  IMSL  Version  10  double  precision  program 
DCSDliR.  The  SFS  technique  was  developed  as  a  separate  generic  computer  program 
subroutine. 

Seven  dilTcrentiable  mathematical  bathymetric  data  models  were  employed  in  deter¬ 
mining  the  accuracy  and  relative  worth  of  both  the  Akima  cubic  spline  and  SI'S  repres¬ 
entations  of  the  ocean  bottom  contour,  including  the  first  and  second-order  derivatives 
of  these  contours.  A  series  of  numerical  enhancements  were  also  applied  to  both  mod¬ 
eling  techniques  to  optimize  their  utility  in  accurately  reconstructing  sampled  ocean 
bottoms. 

The  second  objective  of  this  research  was  the  development  and  evaluation  of  a  re¬ 
flection  angle  algorithm,  derived  from  the  principles  of  the  calculus  as  outlined  by 
Goodman  [Ref  1],  which  enables  the  computation  of  reflection  angles  for  acoustic  rays 
incident  upon  an  ocean  bottom  whose  variation  in  depth  is  modeled  as  a  function  of 
down  range.  This  algorithm  requires  exact  first  derivatives  for  angle  calculations.  Hence 
the  apriori  requirement  for  accurate  ocean  bottom  modeling. 

Both  of  the  developments  outlined  above  were  then  incorporated  into  an  existing 
acoustic  ray  trace  program  and  several  simple  case  studies  investigated. 
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II.  THEORETICAL  ANALYSIS 


A.  SPATIAL  FOURIER  SERIES  REPRESENTATION  OF  OCEAN  BOTTOM 
CONTOURS 

For  any  ocean  region  of  horizontal  range  extent  Z  meters  with  a  one-dimensional, 
range  dependent  bathymetr>’,  the  ocean  bottom  contour  >’^(z)  can  be  modeled  by  the 
following  spatial  Fourier  series  (SFS): 

5 

yeU)  =  I  ^0 1  +  2^  I  Wj  I  cos(27t5/2^z  -  Lw,)  ,  0  <  z  <  Z.  (2.1) 

j=i 

The  magnitude  and  phase  angle  of  the  complex  Fourier  series  coefficient  w„  for  spatial 
harmonic  s,  are  given  by  |  u',|  and  Lw, ,  respectively.  The  average  value  for  depth,  which 
corresponds  to  zero  spatial  frequency  or  the  "DC  component"  of  the  bathymetry  data, 
is  represented  by  |  iVpl .  The  fundamental  spatial  frequency  is  defined  as: 

/  =  (2.2) 

where  N  is  the  total  number  of  bathymetric  data  points  and  Ts^  is  the  uniform  sampling 
period  in  the  z  direction  in  meters  per  sample.  The  highest  spatial  harmonic  is  S  and  is 
governed  by  the  number  of  data  points: 

S  =  —  ,  for  even  N  (2.3) 

and 

5  =  A.tJ-  ^  for  odd  N.  (2.4) 

The  first  and  second-order  derivatives  of  the  SFS  representation  of  discrete  ocean 
bathymetry  are  readily  determined  from  Equation  (2.1): 

s 

3'b'(z)  =  I  I  sin(27rs/2^z  -  Lw^)  ,  0  <  z  <  Z  (2.5) 

j=i 


2 


and 


Xb  ~ 


s 

I 


J=1 


I  I  ^  cos{2nsf2  2  —  L\v^  ,0<z<Z. 


(2.6) 


The  complex  Fourier  series  coelTicients  vv,  are  determined  from  the  spatial  discrete 
Fourier  transform  (DFT)  of  the  ocean  bathymetric  data  as  follows: 


N-\ 


W, 


n=0 


The  term  IV’JJ  is  known  as  the  DFT  weighting  factor  and  is  given  by 


where 


'V=  2S+  1 


(2.7) 


(2.S) 


(2.9) 


when  the  ocean  bottom  sampling  frequency  equals  the  Nyquist  rate. 

The  reconstruction  of  any  one-dimensional  function  or  contour  by  a  Fourier  series 
of  finite  extent  will  be  characterized  by  small  oscilliatory  deviations  from  the  exact 
function  values  in  the  interpolated  region  between  the  discrete  samples.  These  oscil¬ 
lations  are  often  of  little  consequence  due  to  the  small  amplitudes  involved.  However, 
efforts  to  model  a  simple  ocean  bottom  profile  of  constant  non  zero  slope  by  SFS 
methods  generate  significant  amplitude  errors  at  the  function  extremeties.  This  develop¬ 
ment  highlights  a  particular  problem  that  SFS  methods  have  with  apparent  function 
discontinuities  represented  by  the  boundaries  of  the  discrete  data  string.  This  anomoly, 
first  documented  by  Willard  Gibbs  in  1899  and  known  as  "Gibbs'  phenomenon"  [Ref 
2],  can  be  partially  redressed  by  the  inclusion  of  more  data  samples.  This  results  in  the 
compression  of  the  larger  amplitude  oscillations  toward  the  discontinuities.  However,  the 
peak  amplitude  of  these  ripples,  being  independent  of  the  number  of  discrete  data,  re¬ 
mains  constant.  The  problems  associated  with  SFS  methods  and  the  "Gibbs'  phenome¬ 
non"  take  on  greater  significance  when  derivatives  of  the  bottom  contour  are  sought. 
The  ripples  generated  at  the  extremities  by  the  Fourier  series  representation  will  be 
magnified  in  subsequent  calculations  making  reliable  determination  of  function 
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derivatives  difTicult.  Strum  and  Kirk  [Ref.  2]  provide  a  more  detailed  treatment  of 
"Gibbs'  phenomenon". 

Lanezos  [Ref  .3]  argued  that  the  impact  of  this  phenomenon  could  be  minimized  by 
hastening  the  convergence  of  the  Fourier  series  and  proposed  the  following  modified 
complex  Fourier  series  coefficient  as  a  means  of  achieving  this: 


smisnjS) 
snIS 


,  5  7^  0 


(2.10) 


where  w,  is  the  original  coefficient  given  by  Equation  (2.7). 

An  additional  modification  to  Lanezos'  formula  was  suggested  by  Hamming  [Ref 
4]  and  involved  increasing  the  harmonic  S  by  1.  The  final  expression  for  spatial  complex 
Fourier  series  coelficients  employing  the  Lanezos  smoothing  methodology  therefore  be¬ 
comes 


sin(s7r/(S  +  1)) 
57l/(S  +  1) 


Uj  ,  5  #  0. 


(2.11) 


B.  THE  REFLECTION  ANGLE  ALGORITHM 

Acoustic  rays  obey  the  reflection  laws  applied  in  optics.  Equation  (2.12)  relates  the 
angle  of  incidence  C  and  angle  of  reflection  C  for  a  sound  ray  impinging  the  ocean  sur¬ 
face  or  ocean  bottom: 


C  =  C'.  (2.12) 

Consider  a  ray  launched  from  a  sound  source  with  initial  angle  and  propagating 

in  an  isospeed  ocean  with  a  perfectly  flat  sea  surface.  When  the  ocean  bottom  is  flat,  the 
angle  of  incidence  C  will  remain  constant  for  successive  reflections  as  shown  in  Figure 
1.  The  magnitude  of  the  reflection  angle  is  given  by 

C'  =  /Jo.  (2.13) 

When  the  ocean  bottom  shoals  (up  slope)  with  constant  angle  of  slope  y,  the  re¬ 
lationship  given  by  Equation  (2.13)  is  invalidated  and  subsequent  acoustic  ray  reflections 
become  linearly  dependent  upon  the  bottom  slope.  This  observation,  reported  by  Officer 
[Ref  5]  and  illustrated  in  Figure  2,  is  given  by 

C„  =  (/?0-«>-)  ,n=  1,2,3,...  (2.14) 
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where  (/;  =  I)  corresponds  to  the  first  bottom  bounce  with  higher  values  for  n  repres¬ 
enting  subsequent  surface  and  bottom  reflections.  Negative  values  for  s„  indicate  acous¬ 
tic  rays  which  have  been  reflected  back  toward  tlieir  source.  Assigning  a  negative  value 
to  y  in  Equation  (2.1-1)  yields  the  relationship  for  s„  after  an  initial  reflection  from  a 
bottom  contour  that  deepens  (down  slope)  with  constant  angle  of  slope  y.  Figure  3  and 
Equation  (2.15)  refer  to  this  situation; 


C.  =  (/?o +  /')•)  ,n=  1,2,3...  (2.15) 


Figure  1.  Acoustie  Ray  Impinging  a  Flat  Ocean  Butiuin:  fhis  diagram  represents 
an  acoustic  ray  launched  with  initial  angle  of  propagation  /ij,  in  an 
isospeed  ocean  with  a  flat  ocean  surface. 
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Surface 


Figure  2.  Acoustic  Ray  iiupiiiging  a  Constant  IJp-Slope  Ocean  Bottom:  Tliis  dia¬ 
gram  represents  an  acoustic  ray  launched  wiili  initial  angle  of  propa¬ 
gation  /io  in  an  isospeed  ocean  with  a  flat  ocean  surface. 


Figure  3.  Acoustic  Ray  Impinging  a  Constant  Doun-Slope  Ocean  Bottom:  I  his 
diagram  represents  an  acoustic  ray  launched  with  initial  angle  of  propa¬ 
gation  in  an  isospeed  ocean  with  a  Hat  ocean  surface. 
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These  laws  of  reflection  are  equally  valid  in  more  general  ocean  models.  However, 
their  application  is  correspondingly  more  complex.  In  the  case  of  constant  sound  speed 
c,  acoustic  rays  travel  in  straight  lines  within  the  initial  YZ  launch  plane  making  the 
calculation  of  (  relatively  simple.  However,  if  sound  speed  is  made  to  vary  with  depth 
such  that 


c  =  c(y),  (2.16) 

then  the  acoustic  ray  will  remain  in  the  initial  YZ  launch  plane  but  travel  a  curved  path. 
Additionally,  if  the  ocean  bottom  contour  varies  arbitrarily  with  down  range  z  from 
the  sound  source  where 


(2.17) 

then  the  relationships  given  by  Equations  (2.14)  and  (2.15)  are  no  longer  relevant  since 
the  angle  of  slope  y  does  not  remain  constant. 

The  calculation  of  C  in  an  ocean  model  characterized  by  Equations  (2.16)  and  (2.17) 
is  achieved  by  representing  the  direction  of  the  acoustic  ray  at  any  point  along  its  path 
by  the  unit  vector 

/i(y)  =  v{y)j  +  w(>’)  k  (2. IS) 


as  shown  in  Figure  4. 


Figure  4.  Vector  Representation  of  an  Acoustic  Ray;  This  diagram  shows  the 
vector  representation  of  a  sound  ray  travelling  in  an  ocean  medium 
where  sound  speed  is  a  function  of  depth  y. 


7 


The  vector  components  v(y)  and  h  (j-)  are  the  direction  cosines  vv'ith  respect  to  the  Y  and 
Z  axes,  respectively,  describing  the  direction  of  propagation  of  the  ray  and  are  given  by 

v(y)  =  cospiy)  (2.19) 

and 

w(y)  =  sin  /?0).  (2.20) 

In  addition,  the  tangent  plane  and  normal  vector  at  any  point  P{x^,z^)  on  the  ocean 
bottom  must  be  determined. 

The  tangent  plane  at  a  particular  point  on  a  two-dimensional  surface 

given  by 


(2.21) 

can  be  determined  as  follows  (see  Goodman  [Ref.  I]): 

y  -  J'o  =fx{X(i,2f^){x  -  Xq)  -y  f,{xQ,ZQ){z  -  zo)  (2.22) 


where 


(2.23) 


and 


/^(•^0.20)  =  '§^Ax,2)  I  ;c=  ^0-  (2-24) 

Goodman  further  showed  the  unit  vector  normal  to  the  tangent  plane  at  P{x^^\,z^)  to 
be 


A  «x 

=T=n 


(2.25) 


where 


rti  = 


-Xc(^O-?0)'  +J  -f2ix0,zo)k 


(2.26) 
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and 


(2.27) 


Substitution  of  the  function  j.’g  given  by  Equation  (2.17)  into  Equation  (2.22)  yields  the 
following  general  expression  for  the  tangent  plane  at  any  point  PCjo.z©)  on  an  ocean 
bottom  which  is  an  arbitrary  function  of  down  range  z. 


y  -  yo  =  y  -  ^o)- 


(2.28) 


Similarly,  substitution  of  the  function  yg  into  Equation  (2.25)  provides  the  following 
general  relationship  for  the  vector  normal  to  the  bottom  at  any  point  Piy^,!^): 


where 


A 

= 


A  A 

"b1  =  +J  -y’Bi^o)^ 


and 


I  '^Bi  I  =  \/*  +  0''b(’o))^  • 


(2.29) 


(2.30) 


(2.31) 


The  angle  of  slope  y  at  any  point  along  the  one-dimensional  ocean  floor  can  be 
determined  from  elementary  vector  calculus  with  the  use  of  Equation  (2.29)  as  follows: 

y  =  cos“'(v  .«Bj_)  (2.32) 

where  y-  is  the  unit  vector  lying  along  they'  axis.  The  vector  representation  of  an  acoustic 
ray  given  by  Equation  (2.18)  can  be  used  in  defining  the  sound  ray  incident  upon  the 
ocean  bottom  where 


rJi^vJ  +  Wik.  (2.33) 

The  angle  of  arrival  /?„  which  is  measured  from  the  y;  axis  to  «,  and  is  necessary  in  the 
calculation  of  the  angle  of  incidence  [,  is  given  by 

/?.  =  cos'’Cv-«,).  (2.34) 
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Finally,  with  general  expressions  determined  for  y  and  /?„  the  angle  of  incidence  C 
and  the  launch  angle  for  the  reflected  ray  [i,  can  be  evaluated  from  the  following  re¬ 
lationships; 


C  =  /?,--Aj-By  (2.35) 

and 

/?,=  180°-(C  +  A^'bV)  (2.36) 

where 

*  y^fe) 

^^fl  =  TT7Tr-  (2.37) 

\y  I 

Figures  5  and  6  illustrate  the  respective  cases  of  an  acoustic  ray  incident  upon  an  ocean 
bottom  with  negative  and  positive  gradients,  respectively. 
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Figure  5.  Acoustic  Ray  Impinging  an  Arbitrary  Ocean  Buttuin:  I  his  diagram  re¬ 
presents  an  acoustic  ray  incident  upon  an  ocean  bottom  which  is 
shoaling  (up  slope)  at  the  point  of  impact. 
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deepening  (down  slope)  at  the  point  of  impact. 
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III.  OCEAN  BOTTOM  MODELING  OF  DISCRETE  BATHYMETRY 


A.  A  COMPARATIVE  ANALYSIS  OF  AKIMA  CUBIC  SPLINE  AND  SPATIAL 
FOURIER  SERIES  REPRESENTATIONS  OF  BATHYMETRIC  DATA 

The  Akima  Cubic  Spline  (ACS)  algorithm  accurately  models  most  continuous 
functions  from  discrete  data  but  will  frequently  return  inaccurate  first  and  second-order 
derivatives.  Reliable  derivative  calculations  are  crucial  for  reflection  angle  and  ray 
acoustic  computations  and  consequently  an  alternative  modeling  technique  was  needed. 
One  of  the  primary  focuses  of  this  thesis,  therefore,  was  the  investigation  of  Spatial 
Fourier  Series  (SFS)  methods  for  modeling  the  ocean  bottom  from  discrete  bathymetric 
data.  This  particular  approach  possessed  the  added  attraction  in  that  bottom  roughness 
could  be  related  to  spatial  frequency  components.  Furthermore,  derivative  calculations 
could  be  readily  determined  from  the  resultant  series  representation  of  the  contour  ob¬ 
viating  the  requirement  to  call  complex  external  computer  programs. 

Synthetic  bathymetric  profiles  comprising  11  equidistant  data  points  over  a  1000 
meter  range  were  generated  from  seven  exact  mathematical  functions  containing  features 
consistent  with  real  ocean  bottoms.  ACS  and  SFS  representations  of  each  bottom  con¬ 
tour,  including  their  first  and  second-order  derivatives,  were  then  computed.  Percentage 
error  calculations  were  made  with  respect  to  the  exact  function  values  for  each  model 
and  its  associated  derivatives.  When  the  difference  bctw'ccn  the  modeled  and  exact  data 
was  of  magnitude  less  than  10  ‘ ,  the  percentage  error  w’as  assumed  to  be  zero.  Similarly, 
percentage  eriors  greater  than  10‘  were  represented  in  tabular  output  as  a  series  of  as¬ 
terisks. 

The  exact  functions  used  to  generate  the  artificial  ocean  bottom  contours  were; 

•  Flat  bottom  (  Vfi(z)  =  100  m  ) 

•  Constant  downslope  bottom  ( =  220.0  +  0.22  m  ) 

•  Constant  upslope  bottom  ( j;£(z)  =  420.0  —  0.2z  m  ) 

•  One  period  cosine  bottom  ( ^^(z)  =  450.0  —  250.0  cos(0.0027rz)  m  ) 

•  Modified  Witch  of  Agnesi  bottom  (  ygiz)  =  100.0  -  '^'5'j'7/:^5y5  qJ  ) 

•  Half  catenary  bottom  ( >^£(2)  =  3 1 5.0  —  30.0(e^''“‘’  “  +  m  ) 

(z  -  5000.0)' 

•  Parabolic  bottom  (  ^’^(z)  =  50.0  H - ^  ^ -  m  ). 
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Figures  7,  8  and  9  and  Tables  1  and  2  reveal  shortcomings  associated  with  the  SFS 
representation  of  a  flat  ocean  bottom.  While  both  the  ACS  and  SFS  techniques  repli¬ 
cated  precisely  the  original  contour,  the  latter  introduced  small  errors  in  both  the  first 
and  second-order  derivatives. 

Figures  10,  11  and  12  and  Tables  3  and  4  provide  a  comparison  between  the  ACS 
and  SFS  representations  of  an  ocean  bottom  which  deepens  (down  slope)  with  constant 
slope.  Table  3  and  Figure  10  show  that  the  ACS  technique  accurately  modeled  the  ori¬ 
ginal  function  and  its  derivatives.  However,  Figure  10  and  Table  4  highlight  difficulties 
in  the  SFS  representation  of  this  same  contour.  This  method  accurately  reproduced  the 
original  discrete  bathymetric  data  values  but  introduced  oscilliator\’  errors  in  the  in¬ 
terpolated  region  between  them.  These  small  oscillations  around  the  true  function  values 
were  propagated  through  the  subsequent  derivative  calculations  generating  large  per¬ 
centage  errors. 

Figures  13,  14  and  15  and  Tables  5  and  6  provide  an  analogous  comparison  for  an 
ocean  bottom  which  shoals  (up-slope)  at  a  constant  rate.  The  results  observed  for  this 
model  were  identical  to  those  for  the  constant  down  slope  ocean  bottom. 

Figures  16,  17  and  18  and  Tables  7  and  8  pertain  to  the  ACS  and  SFS  represent¬ 
ations  of  the  one  period  cosine  ocean  bottom.  Both  techniques  accurately  reproduced 
the  fundamental  function  shape  as  shown  by  Figure  10  but  the  ACS  approach  had 
trouble  with  both  the  first  and  second-order  derivatives.  Figures  17  and  18  and  Table  7 
refer.  Additionally,  both  the  ACS  and  SFS  reconstructions  of  the  first  and  second-order 
derivatives  were  characterized  by  large  errors  at  the  contour  end  points. 

Figures  19,  20  and  21  and  'fables  9  and  10  provide  a  comparison  between  the  ACS 
and  SFS  representations  of  the  modified  Witch  of  Agnesi  ocean  (  see  [Ref.  6]  )  bottom 
model.  The  SFS  technique  was  incapable  of  reproducing  this  contour  or  its  derivatives 
as  shown  in  the  referenced  Figures  and  Table  10.  Furthermore,  the  oscilliatoiy  errors 
first  observed  with  the  constant  up  slope  and  down  slope  ocean  bottoms  were  again 
present.  The  ACS  approach  yielded  a  far  more  accurate  reproduction  of  the  exact  bot¬ 
tom  function.  However,  both  the  first  and  second-order  derivative  values,  although  bet¬ 
ter  than  those  reproduced  by  SFS  methods,  were  also  unsatisfactory. 

Figures  22,  23  and  24  and  Tables  II  and  12  relate  to  the  ACS  and  SFS  represent¬ 
ations  of  the  half  catenarx'  (  see  [Ref  6] )  ocean  bottom.  The  problems  observed  with  the 
SFS  reconstruction  of  this  contour  were  identical  to  those  observed  for  the  modified 
Witch  of  Agnesi  bottom  as  seen  in  'fable  12.  The  ACS  method  again  provided  a  true 
reproduction  of  the  exact  mathematical  bottom  model.  In  addition,  the  first-order 
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TECHNIQUES  EOR  NUMERICRL  INTERPOLfiTION  OE  YBIM) 


Figure  9.  Flat  Ocean  Buttuin  :  Second-Order  Derivatives. 
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WJTHDin  LRNCZOS  SMOOTHING 


Figure  10.  Down-Slope  Ocean  Bottom  :  Contour  Reconstruction. 
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Figure  11.  Donn-Slupe  Ocean  Botloni  :  First-Order  Derivatives. 
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Figure  13.  Up-Slope  Ocean  Bottom  :  Contour  Reconstruction. 
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Figure  14.  Up-Slope  Ocean  Bottom  :  First-Order  Derivatives. 
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Figure  16.  Cosine  Ocean  Bottom  :  Contour  Reconstruction. 
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Figure  19.  Modified  VVitcli  Of  Agnesi  Ocean  Bottom  :  Contour  Reconstruction. 
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21.  Modified  VV^itch  Of  Agnesi  Ocean  Bottom  :  Second-Order  Derivatives. 
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Table  1.  AKIMA  CUBIC  SPLINE  RESULTS  :  FLAT  OCEAN 
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Table  6.  FOURIER  SERIES  RESULTS  :  CONSTANT  UP-SLOPE 
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derivative  results  were  accurate  with  a  uniformly  small  percentage  error  as  reflected  in 
Table  11.  Second-order  derivative  values  generated  by  the  ACS  approach  were  unsatis¬ 
factory,  although  this  is  not  immediately  evident  from  Figure  24  because  of  the  plot 
scale. 

Figures  25,  26  and  27  and  Tables  13  and  14  give  the  ACS  and  SFS  results  for  the 
parabolic  ocean  bottom  model.  The  ACS  method  matched  both  the  original  function 
and  its  first  and  second-order  derivatives  precisely  as  shown  in  Table  13.  The  SFS  rep¬ 
resentation  was  reasonably  accurate  although  some  oscillations  were  again  in  evidence 
at  the  contour  end  points.  These  oscillations  account  for  the  large  boundary  errors  in  the 
subsequent  first  and  second-order  derivative  computations  as  reflected  in  Table  14  and 
illustrated  by  Figures  26  and  27. 

B.  NUMERICAL  ENHANCEMENT  OF  THE  SPATIAL  FOURIER  SERIES 
TECHNIQUE 

The  results  of  the  previous  section  clearly  showed  that  the  elementary  SFS  approach 
was  unsatisfactory.  Modifications  to  the  fundamental  SFS  technique  were  therefore  in¬ 
vestigated  and  applied  where  appropriate  to  improve  overall  performance. 

The  following  numerical  enhancements  were  studied  for  each  of  the  original  seven 
exact  mathematical  models.  For  purposes  of  brevity,  interim  results  for  the  modified 
Witch  of  Agnesi  contour  only  have  been  included,  as  this  function  reflected  the  greatest 
improvement. 

1.  Lanezos  Sniootliing. 

The  modified  expression  for  the  SFS  coefficients  given  by  Equation  (2.11)  and 
incorporating  Lanezos  smoothing  was  introduced  into  the  SFS  computer  algorithm  and 
the  modeling  process  repeated. 

Figure  28  and  Table  15  show  that  Lanezos  smoothing  was  effective  in  the  sup¬ 
pression  of  the  between  data  point  oscillations  witnessed  previously  in  Table  10  and 
Figure  19.  The  corresponding  reduction  in  the  magnitude  of  error  for  the  first  and  sec¬ 
ond-order  derivatives  can  be  seen  in  Figures  29  and  30.  One  undesirable  effect  of  this 
smoothing  technique  was  the  introduction  of  artificial  curvatures  at  the  end  points  of  the 
bottom  contour.  The  cause  of  this  phenomenon  is  attributed  to  the  apparent  disconti¬ 
nuity  of  the  function  generated  by  the  immediate  commencement  and  cessation  of  data 
at  cither  end. 
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TECHNIQUES  FOR  NUMERICRL  INTERPOLRTION  OF  y(R) 
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Figure  24.  Half  Catenary  Ocean  Buttuni :  Second-Order  Derivatives. 
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Figure  25.  Parabolic  Ocean  Bottom  :  Contour  Reconstruction. 
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Modified  Wilch  Of  Agnesi  Ocean  Bottom  Contour 
Reconstruction:  Fourier  Scries  reconstruction  with  Lanczos  smootli- 


Figure  29 


Modilied 

Derivatives: 


Witch  Of  Agnesi  Ocean  Bottuin  First-Order 
Fourier  Series  reconstruction  with  Lanezos  smoothing. 


Figure  30.  Modil'ied 
Derivatives: 


Witch  Of  Agnesi  Ocean  Bottom  Second-Order 
Fourier  Series  reconstruction  witli  Lanezos  smoothing. 
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Table  11.  AKIMA  CUBIC  SPLINE  RESULTS 
CATENARY  OCEAN  BOITOM 


HALF 


50 


Table  12.  FOURIER  SERIES  RESULTS  :  HALF  CATENARY 
OCEAN  BOnOM 
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Table 


.  AKIMA  CUBIC  SPLINE  RESULTS  :  PARABOLIC 
OCEAN  BO  n  OM 
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Table  N.  FOURIER  SERIES  RESULTS  :  PARABOLIC  OCEAN 
BOITOM 
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Table  15.  FOURIER  SERIES  RESULTS  :  LANCZOS  SMOOTH¬ 
ING. 
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2,  Lanczos  Smoothing  With  Synthetic  Data. 

The  ability  of  Lanczos  smoothing  to  suppress  oscillations  generated  by  the  SFS 
technique  gave  it  great  utility.  However,  the  artificial  curvature  at  the  end  points  of  the 
contour  was  an  undesirable  eflcct  which  needed  to  be  eliminated.  It  was  unlikely  that  the 
problems  associated  with  Lanczos  smoothing  could  be  countered  by  further  modifying 
the  expression  for  the  SFS  coefficients.  Extension  of  the  original  bathymetric  array  by 
adding  equally  spaced  synthetic  data  points  was,  however,  viewed  as  a  possible  solution. 
The  rationale  of  this  procedure  was  to  force  the  undesirable  features  generated  by 
Lanczos  smoothing  beyond  the  actual  region  of  interest. 

Let  each  data  point  be  represented  as  a  function  of  incremented  range  such 

that 


with 


Tfi  =Tb(^-/)  .  »  =  0,1,2,...,(A’ZZ??7'S-  1) 


(3.1) 


z.=  /Az  ,  f  =  0,1,2,...,(A^ZZJ/’7;S  -  1) 


(3.2) 


where  Az  is  the  distance  between  each  data  point  and  NZBPTS  is  the  total  number  of 
original  discrete  data.  Now  let 


/H5  = 


3’/?(^i)  -T/?(^o) 

Az 


(3.3) 


and 


m,  = - - - .  (^.4) 


Synthetic  data  points  were  then  created  by  linear  extrapolation,  that  is. 


j.^b(z_,)  =y’j5(zo)  +  nis  z_i ,  i  =  1,2,...,SFAT///’  (3.5) 


w'here 


(3.6) 


and 


yB{^(SZBPTS-l+i)  ~  y b{'^{NZBPTS-1^  +  ^{NZBPTS-i+t)  >  '  “  ^  ■S 1  A  THP.  (3.7) 
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SYNTH P  is  the  number  of  synthetic  points  to  be  generated.  Figure  31  provides  a 
diagramatic  representation  of  this  methodology. 


Range  of 
Interest 

0 


O  =  Synthetic  Data  Points 


•  =  Original  Bathymetric  Data 

Figure  31.  Generation  Of  Synthetic  Data:  Synthetic  data  points  are  generated 
by  linear  extrapolation. 

Figures  32,  33  and  3“!  and  Table  16  reveal  that  the  addition  of  live  synthetic 
points  to  either  end  of  the  original  data  array  was  successful  in  shifting  the  end  effects 
induced  by  Lanezos  smoothing  outside  the  horizontal  range  of  interest.  This  resulted  in 
a  better  curve  fit  for  the  modified  Witch  of  Agnesi  bottom  and  smaller  errors  in  the 
subsequent  first  and  second-order  derivatives.  One  undesirable  feature  of  the  linear  ex¬ 
trapolation  method  for  simulating  data  was  its  tendency  to  suj>press  the  natural  curva¬ 
ture  of  the  function  being  modeled,  'fhis  phenomenon  is  observed  in  Figure  32. 

Generating  synthetic  data  by  linear  extrapolation  was  a  valid  technique  wheji 
the  nature  of  the  function  being  modeled  was  relatively  featureless.  Ilow'ever,  it  was  not 
suited  to  functions  w'ith  substantial  variation  in  shape.  Consequently,  an  alternative 
method  involving  the  mirroring  of  the  original  data  about  the  end  points  of  the  contour 
was  created.  This  method  is  described  by 

y'B^  -  =yBi^i)  .  '■  =  1,2 . SYMIIP  (3.8) 
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TECHNIOUCS  rOR  NUMCRICRL  3NTERP0LRT10N  OP  YB(M) 
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z 
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Figure  32.  Modified  Witch  Of  Agnesi  Ocean  Bottom  Contour 
Reconstruction:  Fourier  Series  reconstruction  with  Lanczos  smooth¬ 
ing  and  linear  extrapolation. 
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ORic]»#i.  or  drtr  POJ>frs»  jj  nltbct  or  swncric  orm  points  bctori:  fK)  fricRtS 


Figure  33.  Modified  Witch  Of  Agnesi  Ocean  Bottom  First-Order 
Derivatives:  Fourier  Series  reconstruction  with  Lanczos  smoothing 
and  linear  extrapolation. 
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Figure  34.  Modified  VV'itch  Of  Agnesi  Ocean  Bultuin  Second-Order 
Derivatives:  Fourier  Series  reconstruction  with  Lanezos  smooiliing 
and  linear  extrapolation. 
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Table  16.  FOURIER  SERIES  RESUL'IS  :  LINEAR  EXTRAPO- 
LA  I  ED  DATA. 
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and 


d’/j(‘(.VZa/’7'5-l+/>)  ~  yd^(NZOFTS-\-i))  <  ‘  ~  1 >  .V/ ///\  (3.9) 

Figure  35  illustrates  this  mirroring  process. 


Range  of 
Interest 

0 


O  -  Synthetic  Data  Points 


•  s  Original  Bathymetric  Data 

Figure  35.  (Generation  of  .Sjnthetic  Data:  Synthetic  data  is  generated  by  mirror 
imaging  the  original  data  points. 


Figures  36,  37  and  38  and  Table  17  give  the  results  for  the  SI'S  modeling  of  the 
modified  Witch  Of  Agnesi  contour  using  synthetic  data  points  derived  from  the  miiror- 
ing  technique.  Further  improvement  in  the  overall  curve  fit  and  accuracy  of  the  first  and 
second-order  derivatives  was  observed.  However,  some  leveling  of  the  bottom  contour 
curvature  was  still  prevalent.  The  cause  of  this  phenomenon  was  attributed  to  the  ap¬ 
parent  function  discontinuity  arising  from  the  mirroring  process.  This  problem  is  identi¬ 
fied  in  Figure  35  where  mirror  imaging  at  the  right  contour  boundary'  has  created  a  spike 
in  the  data.  In  an  effort  to  compensate  for  this,  the  SFS  technique  has  smoothed  this 
feature,  thereby  suppressing  the  true  function  curvature  at  the  end  points  of  the  range 
of  interest. 

To  resolve  this  problem  of  curvature  suppression,  the  mirror  imaging  teehnique 
given  by  Equations  (3.S)  and  (3.9)  was  altered  as  follows: 
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Figure  36.  Modilied  Wilcli  Of  Agnesi  Ocean  Buttuni  Contour 
Reconslruciion:  Fourier  Series  reconstruction  with  Lanezos  smooth¬ 
ing  and  mirrored  data. 
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Figure  37.  Modified  Witch  Of  Agnesi  Ocean  Bottom  First-Order 
Derivatives:  Fourier  Series  reconsiruciion  with  Lanczos  smootliing 
and  mirrored  data. 
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Figure  38.  Modified  Witch  of  Agnesi  Ocean  Bottom  Second-Order 
Derivatives:  Fourier  Series  reconstruction  with  Lanezos  smoothing 
and  mirrored  data. 
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Table  17.  FOURIER  SERIES  RESULTS  :  MIRROR  IMAGED 
DAIA. 
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.  /•=  1.2 . SYNTllP 


(3.10) 


and 


yBi^{NZBPTS-\+l^  ~ yBi'^(NZBPTS-i)  ^  /—  1,2 . SYNTllP.  (3.11) 

Figure  39  illustrates  this  modilicatioa. 


Range  of 
Interest 


O  =  Synthetic  Data  Points 
9  =  Original  Bathymetric  Data 


Figure  39.  Generation  of  Synthetic  Data:  Synthetic  data  is  generated  from  a 
shifted  inirror  imaging  of  the  original  data. 


Figures  40,  41  and  42  and  Table  18  reflect  only  a  marginal  improvement  in  the 
modified  Witch  of  Agnesi  curve  fit  as  a  consequence  of  using  the  modified  mirror  imag¬ 
ing  procedure.  Additionally,  for  other  exact  mathematical  models  such  as  the  one  period 
cosine  contour,  this  procedure  proved  detrimental,  again  because  of  curvature  sup¬ 
pression. 

An  alternative  means  of  eliminating  function  discontinuities  created  by  mirror 
imaging  was  investigated.  This  problem  was  most  pronounced  when  the  slope  between 
successive  data  points  increased  toward  the  end  points  of  the  contour.  A  new  method 
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Figure  40.  Modil'ied  VVilcli  of  Agnesi  Ocean  Bottuni  Contour 
Reconslruclion:  Fourier  Scries  reconsiruclion  with  Lanczos  smooth¬ 
ing  and  shifted  mirrored  data. 
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Figure  41.  Modified  Witch  Of  Agnesi  Ocean  Bottom  First-Order 
Derivatives:  Fourier  Series  reconstruction  with  Lanczos  smoothing 
and  shifted  mirrored  data. 
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Figure  42.  Modified  Witch  of  Agnesi  Oce.an  Bottom  Second-Order 
Derivatives:  Fourier  Series  reconstruction  with  Lanczos  smoothing 
and  shifted  mirrored  data. 
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Table  18.  FOURIER  SERIES  RESULTS  :  SHIFTED  MIRROR 
IMAGED  DATA. 
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for  creating  synthetic  data  was  therefore  developed  to  take  advantage  of  this  change  in 
contour  slope.  Let 


yai^x)  -yB('o) 

(3.12) 

'^'10  - 

yB(^2)-yB(^i) 

(3.13) 

-  Az 

ysi^OVZBPrS-l))  ~ y b{~(NZBPTS-2^ 

(3.14) 

■  Az 

and 


"'£2  = 


yii  £5-2))  ysi^iNZBPTS-  3 )) 


Aj 


(3.15) 


where  again  Az  is  the  distance  between  each  data  point.  If 

I//J21 1  <  Imiol 


(3.16) 


and,  or 


I  W£2  1  <  1  /»£!  I  , 


(3.17) 


then  synthetic  data  is  generated  at  the  appropriate  end  of  the  profile  as  follows: 


ysi^-i)  =>'£(20)  +  (^’^("ho  -  .  ‘  =  1,2,...,S  YNTIIP  (3.18) 


and,  or 


3’b(‘{A'Z££75-1+/))  ~  yBi^{NZBPTS-\)^  ('^^(^*£1  ~  /A/Jl£))  ,  /  —  \  ,2,...,SY N  'l  II P  (3.19) 


where 


A/?J  =  Wio  ~  ^*21 


(3.20) 


and 


A/n£  =  ?«£,  -  m£2. 


(3.21) 
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A  comparison  of  the  end  points  in  Figures  35  and  43  demonstrate  the  utility  of  this  new 
method  in  preserving  curvature  and  eliminating  apparent  function  discontinuities. 


Range  of 
Interest 


O  =  Synthetic  Data  Points 


•  s  Original  Bathymetric  Data 


Figure  43.  Generation  of  Synthetic  Data:  Synthetic  data  is  generated  by  revers¬ 
ing  the  rate  of  change  of  slope  of  the  original  boundary  data. 


Enhancing  the  original  bathymetric  data  with  synthetic  points  proved  beneficial 
to  the  SFS  modeling  of  the  bottom  contour,  but  no  one  procedure  for  creating  this  ar¬ 
tificial  data  dominated.  Each  of  the  methods  investigated  showed  merit  under  specific 
circumstances.  Consequently,  the  SFS  computer  algorithm  was  modified  to  incorporate 
the  three  methods  outlined  previously.  If 

I  Awl  =  |w,o  -  ^21 1  <  ^  (3.22) 


and,  or 


I  AW£  I  =  I  /«£,  -  W£2  I  <  £  , 


(3.23) 


where  £  is  some  arbitrarily  small  number,  then  the  synthetic  data  was  constructed  by 
linear  extrapolation  as  given  in  Equations  (3.5)  and  (3.7).  However,  if 
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|m,ol<lm2,|  (3.24) 

and,  or 

|m£i  I  <  Im^jl  ,  (3.25) 

then  mirror  imaging  as  given  by  Equations  (3.8)  and  (3.9)  is  used  to  generate  the  data. 
Finally,  if  the  conditions  given  by  Equations  (3.16)  or  (3.17)  hold  then  the  method  em¬ 
ploying  the  reversal  of  the  rate  of  change  of  slope  is  activated  to  simulate  data  at  the 
appropriate  end  point. 

C.  SPATIAL  FOURIER  SERIES  ANALYSIS  OF  THE  AKIMA  CUBIC  SPLINE 
REPRESENTATION  OF  THE  OCEAN  BOTTOM  CONTOUR. 

Results  from  the  SFS  technique  contained  significant  errors  in  the  reconstruction 
of  both  the  original  contour  and  first  and  second-order  derivatives  despite  Lanczos 
smoothing  and  the  optimization  of  methods  for  generating  synthetic  data.  However,  for 
the  more  complex  ocean  bottom  models  such  as  the  modified  Witch  of  Agnesi  and  one 
period  cosine  functions,  SFS  methods  yielded  the  more  accurate  derivative  values. 

The  ACS  method  provided  precise  reproductions  of  the  original  test  functions,  but 
the  subsequent  first  and  second-order  derivative  values  were  reliable  only  for  the  flat, 
constant  slope  and  parabolic  ocean  bottoms.  Efforts  were  therefore  made  to  incorporate 
the  desired  aspects  of  both  methods. 

The  ACS  representation  of  the  ocean  bottom  model  was  discretized  at  ten  times  the 
sampling  frequency  used  to  generate  the  original  bathymetric  data.  Synthetic  data  was 
then  constructed  as  described  in  Section  B  of  Chapter  3,  with  the  exception  that  Az  was 
now  based  upon  this  new  sampling  frequency.  For  example,  where  previously  5  synthetic 
data  points  covered  a  500  m  range,  now  50  synthetic  points  describe  this  same  range.  A 
SFS  representation  of  the  entire  contour  was  then  obtained.  This  approach  was  applied 
to  each  of  the  seven  exact  mathematical  functions. 

Figures  44,  45  and  46  and  Table  19  provide  results  for  the  SFS  representation  of  the 
flat  ocean  bottom  generated  by  the  ACS  method.  The  original  contour  is  replicated 
precisely  and  the  small  errors  observed  in  both  the  first  and  second-order  derivatives 
approach  the  limits  of  machine  precision.  However,  while  these  errors  are  very  small, 
they  are  actually  larger  than  those  generated  by  the  fundamental  SFS  reconstruction  of 
the  flat  bottom  without  Lanczos  smoothing  or  synthetic  data.  Table  2  refers. 
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Figure  -l-l.  Flat  Ocean  Bottom  Contour  Reconstruction:  Spiiiial  Fourier  Series 
representation  of  the  Akima  Cubic  Spline  Hat  bouoin  model. 
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Figure  *16.  Flat  Ocean  Buttoin  Second-Order  Derivatives:  Spatial  Fourier  Series 
representation  of  the  Akima  Cubic  Spline  Hat  bottom  model. 
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Table  19.  FOURIER  SERIES  RESULTS  :  ACS  FLAT  OCEAN 
BOTIOM 
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Figures  47,  48  and  49  and  Table  20  give  the  results  for  the  SFS  representation  of  the 
constant  down-slope  ocean  bottom  generated  by  the  ACS  technique.  The  percentage 
errors  recorded  for  both  the  bottom  contour  and  the  first  and  second-order  derivatives 
are  two  orders  of  magnitude  lower  than  the  initial  results  given  at  Table  4.  This  re¬ 
duction,  although  significant,  is  not  sulTicient  to  eliminate  residual  errors  likely  to  ad¬ 
versely  impact  on  ray  acoustic  computations. 

Figures  50,  51  and  52  and  Table  21  provide  results  for  the  SFS  representation  of  the 
constant  up-slope  ocean  bottom  model.  The  observations  made  for  this  case  are  identical 
to  those  recorded  for  the  constant  down-slope  bottom. 

Figures  53,  54  and  55  and  Table  22  give  results  for  the  SFS  representation  of  the 
one  period  cosine  ocean  bottom.  The  errors  for  both  the  contour  reconstruction  and  its 
first  and  second-order  derivatives  were  larger  than  those  recorded  in  the  earliest  test 
provided  in  Table  8.  Furthermore,  the  sampUng  of  the  ACS  model  for  this  bottom  shape 
not  only  failed  to  assist  SFS  in  more  accurately  reproducing  the  contour,  but  caused  it 
to  replicate  anomolies  in  the  ACS  reconstruction  of  the  first  and  second-order  deriva¬ 
tives. 

Figures  56,  57  and  58  and  Table  23  provide  results  for  the  SFS  representation  of  the 
modified  Witch  of  Agnesi  ocean  bottom.  This  contour  was  used  to  demonstrate  the  de¬ 
velopment  of  successive  numerical  enhancements  to  the  fundamental  SFS  technique,  and 
therefore,  gives  the  best  opportunity  for  a  complete  comparative  analysis.  Significant 
improvement  was  observed  in  the  modeling  of  this  function  with  SFS  methods  employ¬ 
ing  Lanezos  smoothing  and  synthetic  data  created  by  mirror  imaging.  A  comparison  of 
the  results  in  Tables  17  and  23  show  that  further  reductions  in  error  were  achieved  by 
the  SFS  approach  to  modeling  the  ACS  representation  of  the  modified  Witch  of  Agnesi 
bottom.  However  the  recorded  improvements  in  the  first  and  second-order  derivatives 
remain  insufficient  for  the  purposes  of  existing  ray  acoustics  algorithms. 

Figures  59,  60  and  61  and  Table  24  give  results  for  the  SFS  representation  of  the 
half  catenary  ocean  bottom.  Significant  reductions  in  the  percentage  errors  of  the  con¬ 
tour  were  achieved.  For  the  first  and  second-order  derivatives  reductions  in  percentage 
error  of  three  orders  of  magnitude  were  common.  However,  the  residual  errors  remained 
sizeable  and  prohibited  the  use  of  these  results  in  further  computations. 

Figures  62,  63  and  64  and  Table  25  provide  results  for  the  SFS  representation  of  the 
parabolic  ocean  bottom.  The  method  of  using  the  ACS  model  of  the  test  .function  as 
input  to  the  SFS  procedure  again  proved  beneficial  with  significant  reductions  in  the 
magnitude  of  the  percentage  errors  being  observed  for  the  contour.  Reproductions  of  the 
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Figure  47.  Doun-Slope  Ocean  Bottom  Contour  Reconstruction:  Spatial  Fourier 
Scries  reconstruction  of  the  Akima  Cubic  Spline  down-slope  bottom 
model. 
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Figure  48.  Down-Slope  Ocean  Bottom  First-Order  Derivatives:  Spatial  Fourier 
Series  representation  of  the  Akinia  Cubic  Spline  down-slope  bottom 
model. 
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Figure  49.  Doun-SIope  Ocean  Buttuiu  Second-Order  Derivatives:  Spatial  Fourier 
Series  representation  of  the  Akima  Cubic  Spline  down-slope  bottom 
model. 
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TECHNIOUCS  FOR  NUMERICRL  INTERPOLRTION  OF  YB(M) 


Figure  5U.  Up-Slope  Ocean  Bottom  Contour  Reconstruction:  Spatial  I'ouiicr  Se¬ 
ries  representation  of  the  Akima  Cubic  Spline  up-slope  bottom  model. 
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Figure  51.  Up-Slope  Ocean  Bottom  First-Order  Derivatives:  Spatial  Fourier  Sc¬ 
ries  representation  of  the  Akiina  Cubic  Spline  up-slopc  bottom  model. 
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Figure  52.  Up-Slope  Ocean  Buttuni  Second-Order  Derivatives;  Spatial  Fourier 
Series  representation  of  the  Akima  Cubic  Spline  up-slopc  bottom 
model. 
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Figure  53.  One  Period  Cosine  Ocean  Bottom  Contour  Reconstruction:  Spatial 
Fourier  Series  representation  of  the  Akima  Cubic  Spline  cosine  bottom 
model. 
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Figure  5-4.  One  Period  Cosine  Ocean  Bottom  First-Order  Derivati\es:  Spatial 
Fourier  Series  representation  of  the  Akinia  Cubic  Spline  cosine  bottom 
model. 
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Figure  55.  One  Period  Corine  Ocean  Buttuni  Second-Order  Derivalhes:  Spatial 
Fourier  Series  representation  ol'the  Akiina  Cubic  Spline  cosine  bottom 
model. 
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Cubic  Spline  Witch  of  Agnesi  bottom  model. 


88 


Figure  57.  Modified  Witch  of  Agnesi  Ocean  Dotloni  First-Order 
Derivatives:  Spatial  Fourier  Series  representation  of  the  Akiina  Cubic 
Spline  Witch  of  Agnesi  bottom  model. 
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Figure  58.  Modil'icd  Witch  of  Agnesi  Ocean  Uuttuni  Second-Order 
Derivatives:  Spatial  Fourier  Series  representation  of  the  Akima  Cubic 
Spline  Witch  of  Agnesi  bottom  model. 
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Fi  gure  59.  Half  Catenary  Ocean  Bottom  Contour  Reconstruction:  Spatial  l  oui  ier 
Scries  representation  of  the  Akhna  Cubic  Spline  half  catenary  bottom. 
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Figure  60.  Half  Catenary  Ocean  Buttuin  First-Order  Derivatives:  Spatial  Fourier 
Series  representation  of  the  Akiina  Cubic  Spline  half  catenary  bottom. 
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Figure  61.  Half  Catenary  Ocean  Bottom  Second-Order  Derivatives:  Spatial 
Fourier  Series  representation  of  the  Akima  Cubic  Spline  half  catenary 
bottom. 
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Table  20.  FOURIER  SERIES  RESULTS  :  ACS  DOWN-SLOPE 
OCEAN  BOnOM 
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Table  21.  FOURIER  SERIES  RESULTS 
OCEAN  BOITOM 


ACS  UP-SLOPE 
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Table  22.  FOURIER  SERIES  RESULTS  :  ACS  ONE  PERIOD 
COSINE  OCEAN  BOTPOM 
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first  and  second-order  derivatives  were  very  accurate  over  the  major  portion  oftiie  range 
of  interest,  but  were  characterized  by  large  errors  at  the  outer  end  points  of  this  region, 
as  illustrated  in  Figure  64.  This  phenomenon  is  attributed  to  the  reversal  of  slope  method 
for  generating  synthetic  data. 

D.  ALTERNATIVE  INTERIM  MODELING  TECHNIQUES 

Efforts  to  model  arbitrarily  shaped  ocean  bottom  contours  by  SFS  methods  yielded 
unsatisfactory  results  in  every  case  studied  despite  three  numerical  modifications  de¬ 
signed  to  enhance  the  performance  of  this  technique.  ACS  results  for  the  contours  were 
accurate,  but  the  subsequent  derivative  calculations  were  erroneous.  The  need  for  a 
program  which  could  accurately  model  both  the  bottom  and  its  derivatives  remained 
paramount.  Without  this,  the  precise  calculation  of  reflection  angles  for  acoustic  rays 
impinging  an  undulating  sea  floor  could  not  be  realized. 

Motivated  by  the  above  requirements,  attempts  were  made  to  salvage  the  ACS 
technique  by  utilizing  its  accurate  contour  reproduction  and  applying  alternate  numer¬ 
ical  methods  to  obtain  reliable  values  for  the  first  and  second-order  derivatives.  The 
following  central  differencing  algorithms,  extracted  from  [Ref  7],  were  utilized: 

,  ,  ^  ~Tb(^o  +  2/0  +  b(2o  +  h)  —  8j’g(2(,  —  h)  +.Vb(2o  —  2h) 

y  5(^0)  - - TTT -  (2-26) 


and 


-.'•b(~o  +  2/0  +  16y/;(2o  -b  h)  -  30r/;f2o)  -b  1 61^(20  -  h)  -j-gfzp  -  2h) 

m? 


(3.27) 


where  h  is  some  arbitrarily  small  number,  nominally  1/100  of  the  average  spacing 
between  each  original  data  point. 

The  accuracy  of  Equations  (3.26)  and  (3.27)  was  v'erified  in  separate  tests  using  se¬ 
veral  mathematical  functions.  These  relationships  were  then  applied  to  the  ACS  recon¬ 
structions  of  the  seven  exact  mathematical  ocean  bottom  models  listed  in  Section  1  of 
this  chapter.  The  resultant  first  and  second-order  derivatives  were  an  improvement  on 
previous  methods,  but  below  the  precision  expected.  These  poor  results  are  attributed 
to  the  nature  of  the  ACS  program  and  the  manner  in  which  the  interpolating  splines 
between  each  data  point  are  generated. 

Having  abandoned  both  the  SFS  and  ACS  methods  for  modeling  the  ocean  bottom, 
attention  was  focused  upon  a  computer  routine  which  provided  a  divided  differences 
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Figure  62.  Parabolic  Ocean  Bottom  Contour  Reconstruction:  Spatial  Fourier  Sc¬ 
ries  representation  of  the  Akiina  Cubic  Spline  parabolic  bottom. 
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Figure  63.  Parabolic  Ocean  Bottom  First-Order  Derivatives:  Spatial  Fourier  Se¬ 
ries  representation  of  the  Akima  Cubic  Spline  parabolic  bottom. 
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Figure  64.  Parabolic  Ocean  Bottom  Second-Order  Derivatives:  Spatial  Fourier 
Series  rcpreseiuation  of  the  Akiina  Cubic  Spline  parabolic  bottom. 
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Table  25.  l  OURlER  SERIES  RESULTS  :  ACS  PARABOLIC 
OCEAN  BorroM 
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interpolating  pohnoniial  fit  to  unordered,  randomly  spaced  data.  This  program  ,  written 
by  Gerald  and  Wheatley  [Ref.  7],  was  tested  on  the  seven  exact  mathematical  bottom 
models.  Accurate  curve  fits  were  achieved;  however,  more  care  was  required  in  sampling 
each  contour  to  avoid  artificial  curvatures  in  the  fit  which  are  often  characteristic  of  this 
numerical  technique.  The  subsequent  first  and  second-order  derivatives  generated  by 
applying  Equations  (3.26)  and  (3.27)  to  the  polynomial  fit  were  precise.  Therefore,  these 
equations  and  the  polynomial  coefficients  program  formed  the  basis  of  an  interim  mod¬ 
eling  technique  which  would  enable  verification  of  the  reflection  angle  algorithm  for  an 
arbitrarily  shaped  ocean  bottom. 
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IV.  ACOUSTIC  RAY  REFLECTIONS  FROM  ARBITRARILY  SHAPED 

BOTTOM  CONTOURS 


The  reflection  angle  algorithms  given  by  Equations  (2.13),  (2.14),  (2.15)  and  (2.35) 
were  incorporated  into  an  existing  acoustic  ray  trace  program.  This  Fortran  program 
had  previously  been  limited  to  the  flat  bottom  case.  Seven  exact  mathematical  functions 
where  then  used  to  model  ocean  bathymetry.  One  acoustic  ray  was  launched  for  each 
bottom  contour  and  the  computer  generated  angles  of  reflection  and  bottom  slope 
compared  with  the  corresponding  theoretical  values.  An  isospeed  ocean  was  assumed  in 
each  case  in  order  to  interpret  the  results  easily. 

The  seven  exact  mathematical  functions  employed  were; 

•  Flat  Bottom  (y^(2)  =  1000  m) 

•  Constant  down-slope  bottom  (y^(z)  =  600.0  +  0.042  m) 

•  Constant  up-slope  bottom  (^^(z)  =  1000.0  —  0.04z  m) 

•  One  period  cosine  bottom  (ysi^)  =  250.0  —  50.0  cos(0.0002;rz)  m) 

•  Modified  Witch  of  Agnesi  bottom  0'i<(2)  =  200.0  -  q  o00U^^^+^2500"o' 

•  Half  catenary  bottom  Cy2((2)  =  315.0  —  30.0(e'/‘’““''*’  +  m) 

(z  -  5000.0)2 

•  Parabohc  bottom  (>'s(z)  =  250.0  -I - 7Jq~oo"5 — 

Figures  65  through  7S  illustrate  the  ocean  bottom  models  and  the  propagated 
acoustic  ray  for  each  of  the  test  cases.  A  comparison  between  theory  and  the  computer 
generated  results  for  each  model  are  given  in  Tables  26  through  32  inclusively.  For  each 
table,  the  terms  y°  and  C°  represent  the  exact  values  for  the  angle  of  bottom  slope  and 
acoustic  ray  reflection  angle,  respectively.  The  terms  y°  and  C°  represent  the  estimates 
generated  by  the  ray  trace  program  and  reflection  angle  algorithm  for  these  same  pa¬ 
rameters.  The  bottom  and  surface  bounces  are  denoted  by  BB  and  SB,  respectively. 

These  results  were  uniformly  accurate  with  almost  all  of  the  recorded  errors  for  the 

A 

computed  reflection  angles,  C",  below  0.1%.  Errors  in  the  estimated  bottom  slope,  y°, 
were  also  reasonably  small  thus  verifying  the  central  diflerencing  methods  for  determin¬ 
ing  the  first-order  derivatives.  The  parabolic  ocean  bottom  model  yielded  the  worst  re¬ 
sults  with  errors  in  the  reflection  angle,  C”,  averaging  0.0S%.  The  modified  Witch  of 
Agnesi  and  half  catenary  models  generated  the  largest  errors  in  angle  of  slope,  y°,  but 
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the  magnitudes  of  these  angles  were  sufliciently  small  such  that  little  impact  on  the 
computed  angles  of  reflection  resulted. 

Arithmetic  confirmation  of  the  reflection  angle  algorithm  is  not  readily  determined 
for  a  ray  launched  in  an  ocean  medium  with  speed  of  sound  varying  as  a  function  of 
depth  and  for  an  arbitrarily  shaped  ocean  bottom.  However,  two  such  cases  were  inves¬ 
tigated  as  a  means  of  visual  verification  of  the  effects  of  a  ray  incident  upon  an  arbi¬ 
trarily  shaped  ocean  bottom.  The  one  period  cosine  and  modil'ied  Witch  of  Agnesi 
models  were  utilized  after  slight  modification  to  limit  the  range  of  depth.  Figures  79  and 
80  refer.  The  sound-speed  profile  is  given  by  Table  33. 


RECONSTRUCTED  OCEAN  BOTTOM  BATHTMETRT 


-2.0  -1.0  0.0  1.0  2.0  3.0  1.0 
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Figure  65.  Flat  Ocean  Buttoin:  This  profile  was  used  to  lest  the  reflection  angle 
algorithm  for  an  isospeed  ocean. 
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in  an  isospced  ocean  with  a  flat  bottom. 
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Figure  67.  Constant  Down-Slope  Bottom:  This  profile  was  used  to  test  the  re- 
fiection  angle  algorithm  for  an  isospeed  ocean. 
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Figure  68.  Acoustic  Ray  Impinging  a  Constant  Down-Slope  Bottom:  Tiiis  figure 
sliows  an  acoustic  ray  launched  at  depth  15  ni  with  initial  angle 
j^a  =  20°  and  propagating  in  an  isospeed  ocean  with  a  constant  down- 
slope  bottom. 
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Table  26.  REFLECTION  ANGLES  :  FLAT  DO  R  OM  WITH  /?„  =  45’ 


BOUNCE 

Range 

(kjn) 

ly°l 

ly°l 

%  Error 

A 

c° 

%  Error 

I  BB 

0.9S5 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

1  SB 

I.9S5 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

2  BB 

2.9S5 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

2  SB 

3.985 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

3  BB 

4.985 

0.000 

o.ooo 

0.00 

45.000 

45.000 

0.00 

3  SB 

5.985 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

4  BB 

6.985 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

4  SB 

7.985 

o.ooo 

0.000 

0.00 

45.000 

45.000 

0.00 

5  BB 

8.985 

0.000 

0.000 

0.00 

45.000 

45.000 

0.00 

5  SB 

9.985 

o.ooo 

0.000 

0.00 

45.000 

45.000 

0.00 

Table  27. 


KEFLEC 


ION  ANGLES  :  DOWN-SLOPE  BO  R  OM  WITH  /?„  =  20 


O 


BOUNCE 

Range 

(kill) 

ly°l 

ly°l 

%  Error 

c 

A 

c° 

%  Error 

1  BB 

0.216 

2.291 

2.291 

0.00 

22.291 

22.291 

0.00 

1  SB 

0.494 

0.000 

0.000 

0.00 

24.581 

24.581 

0.00 

2  BB 

0.783 

2.291 

2.291 

0.00 

26.872 

26.872 

0.00 

2  SB 

1.136 

0.000 

0.000 

0.00 

29.162 

29.162 

0.00 

3  BB 

1.504 

2.291 

2.291 

0.00 

31.453 

31.453 

0.00 

3  SB 

1.945 

0.000 

0.000 

0.00 

33.744 

33.744 

0.00 

4  BB 

2.410 

2.291 

2.291 

0.00 

36.0.34 

33.034 

0.00 

4  SB 

2.961 

O.OOO 

0.000 

0.00 

38.325 

38.325 

0.00 

5  BB 

3.547 

2.291 

2.291 

o.oo 

40.615 

40.615 

0.00 

5  SB 

4.237 

0.000 

0.000 

0.00 

42.906 

42.906 

0.00 

6  BB 

4.979 

2.291 

2.291 

0.00 

45.197 

45.197 

0.00 

6  SB 

5.851 

O.OOO 

0.000 

0.00 

47.487 

47.487 

0.00 

7  BB 

6.803 

2.291 

2.291 

0.00 

49.778 

49.778 

0.00 

7  SB 

7.922 

0.000 

0.000 

0.00 

52.069 

52.069 

0.00 

8  BB 

9.162 

2.291 

2.291 

0.00 

54.359 

54.395 

0.00 
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Figure  69.  Constant  Up-Slope  Bottom:  This  profile  was  used  to  test  the  reflection 
angle  algorithm  for  an  isospeed  ocean. 
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Pa  =  55°  and  propagating  in  an  isospced  ocean  with  a  constant  up-slope 
bottom. 
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Table  28.  REFLECTION  ANGLES  :  UP-SLOPE  BOTT OM  WITH  /?„  =  55“ 


BOUNCE 

Range 

(kin) 

ly^l 

ly^l 

%  Error 

A 

c 

%  Error 

1  BB 

0.947 

2.291 

2.291 

0.00 

52.709 

52.709 

0.00 

1  SB 

2.476 

0.000 

0.000 

0.00 

50.419 

50.419 

0.00 

2  BB 

3.515 

2.291 

2.291 

0.00 

48.128 

48.128 

0.00 

2  SB 

4.400 

0.000 

0.000 

0.00 

45.838 

45.838 

0.00 

3  BB 

5.215 

2.291 

2.291 

0.00 

43.547 

43.547 

0.00 

3  SB 

5.909 

0.000 

0.000 

0.00 

41.256 

41.256 

0.00 

4  BB 

6.557 

2.291 

2.291 

0.00 

38.966 

38.966 

0.00 

4  SB 

7.106 

0.000 

0.000 

0.00 

36.675 

36.675 

0.00 

5  BB 

7.624 

2.291 

2.291 

0.00 

34.385 

34.385 

0.00 

5  SB 

8.059 

0.000 

0.000 

0.00 

32.094 

32.094 

0.00 

6  BB 

8.474 

2.291 

2.291 

0.00 

29.803 

29.803 

0.00 

6  SB 

8.818 

0.000 

0.000 

0.00 

27.513 

27.513 

0.00 

7  BB 

9.149 

2.291 

2.291 

0.00 

25.222 

25.222 

0.00 

7  SB 

9.417 

0.000 

0.000 

0.00 

22.931 

22.931 

0.00 

8  BB 

9,676 

2.291 

2.291 

0.00 

20.641 

20.641 

0.00 

S  SB 

9.879 

0.000 

0.000 

0.00 

18.350 

18.350 

0.00 
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propagating  in  an  isospeed  ocean. 
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Table  29.  REFLECTION  ANGLES  :  COSINE  BOTTOM  WITH  =  30“ 


BOUNCE 

Range 

(kjn) 

lv°l 

ly°l 

%  Error 

c° 

O 

< 

%  Error 

1  BB 

0.107 

0.121 

0.124 

2.48 

30.121 

30.124 

0.010 

1  SB 

0.224 

0.000 

0.000 

0.00 

30.242 

30.248 

-0.019 

2  BB 

0.341 

0.383 

0.381 

-0.52 

30.624 

30.628 

0.013 

2  SB 

0.462 

0.000 

0.000 

0.00 

31.007 

31.009 

0.006 

3  BB 

0.584 

0.646 

0.644 

-0.31 

31.652 

31.653 

0.003 

3  SB 

0.713 

0.000 

0.000 

0.00 

32.298 

32.297 

-0.003 

4  BB 

0.843 

0.910 

0.910 

0.00 

33.208 

33.207 

-0.003 

4  SB 

0.983 

0.000 

0.000 

0.00 

34.118 

34.118 

0.000 

5  BB 

1.127 

1.171 

1.173 

0.17 

35.288 

35.290 

0.006 

5  SB 

1.284 

0.000 

0.000 

0.00 

36.456 

36.463 

0.019 

6  BB 

1.446 

1.419 

1.420 

0.07 

37.878 

37.883 

0.013 

6  SB 

1.625 

0.000 

0.000 

0.00 

39.297 

39.302 

0.013 

7  BB 

1.813 

1.634 

1.632 

-1.22 

40.931 

40.935 

0.009 

7  SB 

2.023 

0.000 

0.000 

0.00 

42.566 

42.567 

0.002 

8  BB 

2.245 

1.777 

1.775 

-0.11 

44.342 

44.342 

0.000 

8  SB 

2.497 

0.000 

0.000 

0.00 

46.118 

46.117 

-0.002 

9  BB 

2.766 

1.774 

1.776 

0.11 

47.893 

47.893 

0.000 

9  SB 

3.070 

0.000 

0.000 

0.00 

49.667 

49.668 

0.002 

10  BB 

3.396 

1.522 

1.522 

0.00 

51.189 

51.190 

0.002 

10  SB 

3.759 

0.000 

0.000 

0.00 

52.711 

52.713 

0.004 

11  BB 

4.144 

0.922 

0.922 

0.00 

53.633 

53.634 

0.002 

11  SB 

4.555 

0.000 

0.000 

0.00 

54.555 

54.556 

0.002 

12  BB 

4.977 

0.026 

0.026 

0.00 

54.581 

54.582 

0.002 

12  SB 

5.399 

0.000 

0.000 

0.00 

54.608 

54.608 

0.000 

13  BB 

5.812 

0.879 

0.879 

0.00 

53.728 

53.729 

0.002 

13  SB 

6.200 

0.000 

0.000 

0.00 

52.849 

52.850 

0.002 

14  BB 

6.566 

1.499 

1.498 

-0.07 

51.350 

51.352 

0.004 

14  SB 

6.896 

0.000 

0.000 

0.00 

49.851 

49.854 

0.006 
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Figure  73.  Modified  Witch  Of  Agnesi  Ocean  Bottom:  'I'liis  profile  was  used  to 
test  the  rellection  angle  algorilhin  for  an  isospeed  ocean. 
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/?o  =  70°  and  propagating  in  an  isospced  ocean  with  a  Witch  of  Agnesi 
bottom. 
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Table  30.  REFLECTION  ANGLES  :  WITCH  OF  AGNESI  BOTTOM  WITH 


Po  =  VO’ 


BOUNCE 

Range 

(km) 

ly“l 

ly“l 

%  Error 

O 

< 

%  Error 

1  BB 

0.372 

0.084 

0.084 

0.00 

70.084 

70.084 

0.000 

1  SB 

0.788 

0.000 

0.000 

0.00 

70.168 

70.168 

0.000 

2  BB 

1.212 

0.248 

0.247 

-0.40 

70.416 

70.416 

0.000 

2  SB 

1.647 

0.000 

0.000 

0.00 

70.644 

70.663 

0.027 

3  BB 

2.096 

0.348 

0.348 

0.00 

71.012 

71.011 

-0.001 

3  SB 

2.563 

0.000 

0.000 

0.00 

71.359 

71.359 

0.000 

4  BB 

3.048 

0.371 

0.372 

0.27 

71.730 

71.731 

0.001 

4  SB 

3.554 

0.000 

0.000 

0.00 

72.102 

72.102 

0.000 

5  BB 

4.080 

0.337 

0.336 

-0.30 

72.439 

72.438 

-0.001 

5  SB 

4.629 

0.000 

0.000 

0.00 

72.776 

72.774 

-0.003 

6  BB 

5.196 

0.275 

0.276 

0.36 

73.051 

73.049 

-0.003 

6  SB 

5.784 

0.000 

0.000 

0.00 

73.326 

73.325 

-0.001 

7  BB 

6.388 

0.211 

0.213 

0.95 

73.537 

73.538 

0.001 

7  SB 

7.009 

0.000 

0.000 

0.00 

73.749 

73.752 

0.004 

8  BB 

7.644 

0.157 

0.155 

-1.27 

73.906 

73.907 

0.001 

8  SB 

8.292 

0.000 

0.000 

0.00 

74.063 

74.062 

-0.001 

9  BB 

8.950 

0.116 

0.1 15 

-0.86 

74.179 

74.177 

-0.003 

9  SB 

9.619 

0.000 

0.000 

0.00 

74.295 

74.291 

-0.005 
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Figure  75.  Half  Catenary  Ocean  Buttuni:  This  profile  was  used  to  lest  tlie  rc- 
fiection  angle  algorithm  for  an  isospeed  ocean. 
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Table  31.  REFLECTION  ANGLES  :  HALF  CATENARY  BOTTOM  WITH 


k  =  60° 


BOUNCE 

Range 

(kin) 

ly°l 

ly»l 

%  Error 

A 

%  Error 

1  BB 

0.416 

0.039 

0.035 

-10.25 

59.960 

59.965 

0.008 

1  SB 

0.856 

0.000 

0.000 

0.00 

59.921 

59.931 

0.017 

2  BB 

1.294 

0.125 

0.125 

0.00 

59.800 

59.805 

0.008 

2  SB 

1.728 

0.000 

0.000 

C.OO 

59.671 

59.680 

0.015 

3  BB 

2.157 

0.210 

0.214 

1.90 

59.461 

59.466 

0.008 

3  SB 

2.579 

0.000 

0.000 

0.00 

59.251 

59.253 

0.003 

4  BB 

2.995 

0.298 

0.300 

0.67 

58.953 

58.953 

0.000 

4  SB 

3.401 

0.000 

0.000 

0.00 

58.654 

58.653 

-0.002 

5  BB 

3.799 

0.388 

0.388 

0.00 

58.267 

58.265 

-0.003 

5  SB 

4.186 

0.000 

0.000 

0.00 

57.879 

57.877 

-0.003 

6  BB 

4.563 

0.479 

0.478 

-0.21 

57.400 

57.399 

-0.002 

6  SB 

4.926 

0.000 

0.000 

0.00 

56.922 

56.921 

-0.002 

7  BB 

5.280 

0.572 

0.571 

-0.17 

56.350 

56.350 

0.000 

7  SB 

5.618 

0.000 

0.000 

0.00 

56.922 

56.921 

-0.002 

S  BB 

5.946 

0.665 

0.665 

0.00 

55.113 

55.114 

0.002 

8  SB 

6.258 

0.000 

0.000 

0.00 

54.447 

54.448 

0.002 

9  BB 

6.560 

0.756 

0.759 

0.39 

53.688 

53.689 

0.002 

9  SB 

6.845 

0.000 

0.000 

0.00 

52.929 

52.930 

0.002 

10  BB 

7.120 

0.851 

0.851 

0.00 

52.078 

52.078 

0.000 

10  SB 

7.378 

0.000 

0.000 

0.00 

51.227 

51.227 

0.000 

11  BB 

7.627 

0.941 

0.941 

0.00 

50.287 

50.286 

-0.002 

11  SB 

7.859 

0.000 

0.000 

0.00 

49.346 

49.345 

-0.002 

12  BB 

8.083 

1.027 

1.027 

0.00 

48.319 

41.318 

-0.002 

12  SB 

8.290 

0.000 

0.000 

0.00 

47.291 

47.291 

0.000 

13  BB 

8.490 

1.110 

1.109 

0.09 

46.182 

46.182 

0.000 

13  SB 

8.675 

0.000 

0.000 

0.00 

45.072 

45.073 

0.002 

14  BB 

8.852 

1.187 

1.187 

0.00 

43.885 

43.887 

0.005 

14  SB 

9.015 

0.000 

0.000 

0.00 

42.698 

42.700 

0.005 
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Table  32.  REFLECTION  ANCLES  :  PARABOLIC  BOTTOM  WITH  /?„  =  70“ 


BOUNCE 

Range 

(kin) 

ly°l 

ly°l 

%  Error 

A 

c° 

%  Error 

1  BB 

1.180 

5.816 

5.798 

0.31 

64.184 

64.202 

-0.028 

1  SB 

1.903 

0.000 

0.000 

0.00 

58.368 

58.405 

0.063 

2  BB 

2.451 

3.889 

3.888 

-0.03 

54.479 

54.517 

0.069 

2  SB 

2.861 

0.000 

0.000 

0.00 

50.590 

50.628 

0.075 

3  BB 

3.217 

2.722 

2.721 

-0.04 

47.868 

47.907 

0.081 

3  SB 

3.511 

0.000 

0.000 

0.00 

45.145 

45.186 

0.091 

4  BB 

3.783 

1.859 

1.860 

0.05 

43.287 

43.326 

0.09 

4  SB 

4.021 

0.000 

0.000 

0.00 

41.428 

41.467 

0.094 

5  BB 

4.249 

1.148 

1.147 

-0.09 

40.280 

40.319 

0.097 

5  SB 

4.459 

0.000 

0.000 

0.00 

39.133 

39.172 

0.099 

6  BB 

4.664 

0.514 

0.513 

-0.19 

38.619 

38.658 

0.101 

6  SB 

4.861 

0.000 

0.000 

0.00 

38.104 

38.145 

0.108 

7  BB 

5.057 

0.088 

0.088 

0.00 

38.192 

38.232 

0.105 

7  SB 

5.255 

0.000 

0.000 

0.00 

38.280 

38.320 

0.104 

8  BB 

5.454 

0.695 

0.694 

-0.14 

38.975 

39.014 

0.100 

8  SB 

5.665 

0.000 

0.000 

0.00 

39.670 

39.709 

0.098 

9  BB 

5.881 

1.346 

1.346 

0.00 

41.016 

41.055 

0.095 

9  SB 

6.118 

0.000 

0.000 

0.00 

42.362 

42.401 

0.092 

10  BB 

6.370 

2.092 

2.093 

0.05 

44.454 

44.494 

0.090 

10  SB 

6,661 

0.000 

0.000 

0.00 

46.546 

46.587 

0.088 

11  BB 

6.980 

3.022 

3.020 

-0.07 

49.568 

49.607 

0.078 

11  SB 

7.376 

0.000 

0.000 

0.00 

52.591 

52.627 

0.068 

12  BB 

7.844 

4.337 

4.343 

0.14 

56.928 

56.970 

0.073 

12  SB 

8.498 

0.000 

0.000 

0.00 

61.265 

61.313 

0.078 

13  BB 

9.433 

6.742 

6.669 

-1.08 

68.007 

67.982 

-0.037 
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Figure  19.  Acoustic  Ray  Impinging  A  Cosine  Bottom:  This  ligure  shows  an 
acoustic  ray  launched  at  depth  5  m  with  initial  angle  =  85°  and 
propagating  in  an  ocean  medium  where  sound  speed  varies  as  a  func¬ 
tion  of  depth. 
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/?„  =  85“  and  propagating  in  an  ocean  medium  where  speed  of  sound  is 
a  function  of  depth. 
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Table  33.  SOUND  SPEED  PROFILE 


Depth  (in) 

Sound  Speed  (m/s) 

0.00 

1540.00 

25.00 

1524.00 

30.00 

1522.00 

36.00 

1524.00 

40.00 

1529.00 

45.00 

1532.00 

50.00 

1531.00 

60.00 

1.525.75 

100.00 

1514.75 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  IMSL  Version  10  Akima  Cubic  Spline  (ACS)  computer  computer  program  had 
proven  useful  in  modeling  ocean  bottom  contours,  but  was  unable  to  provide  accurate 
first  and  second-order  derivatives  for  the  one  period  cosine  and  modified  Witch  of 
Agnesi  ocean  bottom  models.  The  application  of  Spatial  Fourier  Series  (SFS)  methods 
as  an  alternative  means  of  modeling  the  one-dimensional  ocean  floor  was  considered  a 
viable  solution  to  this  problem.  However,  despite  several  numerical  enhancements  to  the 
fundamental  technique,  including  Lanczos  smoothing  and  artificial  data  simulation  to 
eliminate  end  point  anomalies.  Spatial  Fourier  Series  failed  to  achieve  the  accuracy  re¬ 
quired  for  subsequent  acoustic  ray  and  reflection  angle  computations. 

Despite  the  aforementioned  difficulties.  Spatial  Fourier  Series  methods  may  yet 
prove  to  be  a  viable  numerical  tool  in  modeling  ocean  bathymetries,  but  its  immediate 
utility,  based  upon  the  analysis  conducted,  is  doubtful. 

The  divided  differences  polynomial  curve  fitting  Fortran  routine  yielded  results 
which  were  equal  to,  or  better  than,  both  the  ACS  and  SFS  modeling  techniques  that  it 
replaced.  Similarly,  the  central  differencing  method  for  determining  the  first  and 
second-order  derivatives  of  the  ocean  bottom  provided  the  most  accurate  estimates  when 
used  in  conjunction  with  this  alternate  approach  to  bottom  contour  modeling. 
Notwithstanding  the  foregoing  observations,  this  polynomial  routine  also  had  deficien¬ 
cies,  the  most  notable  being  the  introduction  of  artificial  oscillations  generated  following 
an  imprudent  choice  of  discrete  bathymetric  data  to  describe  the  contour. 

The  reflection  angle  algorithm  developed  in  Chapter  2  for  an  arbitrarily  shaped 
one-dimensional,  ocean  bottom  were  validated  by  the  computer  simulations  given  in 
Chapter  3. 

One  of  the  realizations  of  this  research  was  that  no  one  numerical  modeling  method 
studied  was  definitive.  Further  investigation  of  bottom  contour  modeling  techniques, 
such  as  the  method  of  least  squares,  is  warranted.  This  work  represents  precursor}'  ef¬ 
forts  in  the  development  of  a  two-dimensional  ocean  bottom  model  which  would  enable 
more  advanced  study  of  ray  acoustic  theor\'  in  an  ocean  medium  where  sound-speed 
varies  as  a  function  of  all  three  spatial  coordinates. 
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